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Abstract 

The suppression of order parameter fluctuations at the boundaries 
causes the ultrasonic attenuation near the superfluid transition to be 
lowered below the bulk value . We calculate explicitly the first devi- 
ation from the bulk value for temperatures above the lambda point 
. This deviation is significantly larger than for static quantities like 
the thermodynamic specific heat or other transport properties like the 
thermal conductivity . This makes ultrasonics a very effective probe 
for finite size effects . 



PACS number(s): 64.60.Ht 



Critical phenomena in confined geometry has been attracting a fair 
amount of attention of late because of the progress on the experimental 
front [2 — 9] which is making it possible to check the predictions of finite size 
effects (FSE). A fair amount of this experimental effort has gone in studying 
the specific heat near the superfluid transition . With the bulk specific heat 
quite well understood and the existence of a sharp phase transition (apart 
from gravity rounding , which too can be removed by doing experiments 
in space) established , efforts have been made to study the FSE . It is ex- 
pected that the FSE will round out the transition and hence the divergence 
at T = T\ will be removed . The specific heat will be finite and the finite 
value will be a function of the confining length . We will keep in mind one 
of the favoured experimental geometries, where one takes two parallel plates 
separated by a distance L , much smaller than the linear dimensions of the 
plates . For L ^> £ , the correlation length at a given temperature, usual 
thermodynamic result follows . It is when L < £ , that FSE dominate . 
Finite size scaling suggests the existence of a scaling function , function of 
£/L- in terms of which the theory can be cast . The specific heat C(t,L) 
in finite geometry has the form C{t,L) ~ t~ a g(t~ u / L) + Constant , where 
£ ~ t~ v and t — (T — T\)/T\ , T\ being the transition temperature . The 
function g has been calculated by various authors [10 — 12] . In what follows 
we propose a method of checking for FSE by studying a related dynamic 
property . This is the study of ultrasonic attenuation (UA) near T\ at high 
frequencies . In fact, it is our contention that UA is one of the best ways 
of checking for FSE since the single surface effect alone can produce effects 
greater than 10% . The critical fluctuations relax according to £~ 2 , where z is 
the dynamic scaling exponent . For frequencies uj such that uj ^> To (occurs 
if one is close to the critical point) , the attenuation is independent of the 
correlation length. For a finite size system , we will show below that in this 
limit , the attenuation is determined by uj and L alone . We provide explicit 
answers for frequencies uj which are much smaller than a cutoff frequency ujq 
(of the order of a few GHz) and for plate separation L > (2r /c^) 2 for a given 
frequency uj . Our prediction for the attenuation per wavelength as a func- 
tion of uj for the plate separation of 2110A is shown in Fig(l) . It should be 
possible to check the prediction experimentally . In fact , this should be the 
simplest way of checking for FSE since the effect is quite pronounced (about 
18% at L — 2110A and uj = lOMhz) for the available confining lengths as 
shown in Fig(l) . This occurs because the imaginary part of the specific heat 
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determines the UA and which is much smaller than the real part , but as we 
shall see below , both are equally affected by the FSE . Consequently, 
the relative effect is much larger for the imaginary part and this will show 
up in the UA . 

The basis of our calculation is once more the Pippard Buckingham Fair- 
bank (PBF) relation [13 — 14] which gives a successful account [15 — 17] of 
the critical ultrasonics in the situation where L ^> £ . The PBF relation is 
obtained from general considerations of entropy clamping and yields for the 
sound velocity u(T,lu) 

u(T, u) = u (T ) + Ul Co/C P (T, cu) (1) 

where Uq(T ) is the sound speed at the transition point (T is the bulk T\ 
for the infinite system , but is a L-dependent temperature for the finite size 
system) , u\ and C are constants and Cp(T, u) is the specific heat at finite 
frequency . 

For the bulk case , C p (T,uj = 0) diverges at T = T A and C p (T,oj) is a 
homogeneous function of cu and £ . If the characteristic relaxation rate is 
To£~ z then the scaling form of C P ulk is 

CF k (T,u)=e*' v f{=^) (2) 

The exponent a is very close to zero for the superfluid transition in in A He 
and for many practical purposes , it is possible to write 

C b P ulk (T,u;) = C[HA0 + f(-^)] (3) 

The function /(cj/(r £ -2 )) reduces to a constant for u = and tends to 
- ln(cj/r ) 1/z £ for uo > T ^ z . A one loop calculation of the scaling function 
f(fl) where Q = r ^- z , was carried out and led to a successful scaling theory 
of the attenuation in the bulk 4 He near T\ [15 — 17] . 

We now need to discuss the effect of a confining geometry . At zero 
frequency , the specific heat is blunted due to the FSE and the usually 
divergent specific heat remains finite . The single loop calculation of the 
scaling function g{^/L) discussed before gives a very reasonable account of 
the recent specific heat data by Mehta and Gasparini [2] . One of the most 
important feature of the scaling function is the low £/L limit (experimentally 
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most easily accessible) is the first departure from the thermodynamic limit 
- the magnitude of this departure AC has to be proportional to the surface 
(A) to volume (V) ratio and hence from purely dimensional arguments , the 
correction can be written as 

AC = C(£, L) - Coo(0 = -aCA^ (4) 

where a is a number of 0(1) , that can be obtained from the function g(£/L) 
, and C is the dimensional constant defined in eqn(3) . The value of a as 
inferred from Scmolke et al. [11] is 1.4. The agreement of this departure 
with the measured departure of Mehta and Gasparini is impressive . 

In our present concern we need the three variable function C(C,,L,uj) 
, whose two limits C(£,u;) and C(£,L) are already well known . We will 
characterize C(£,L, uj) by its first departure from the infinite volume limit 
C(£, u;) and write the generalization of eqn(4) as 

AC(£, L, u ) = C(£, L, u) - C(£, u) = -a(e, v)C(t)A/V (5) 

where a(£, a;) is a scaling function , whose zero frequency limit is a£ (see 
Eqn(4)) and whose general form will be presented below . As soon as we start 
discussing the scaling function for C(£,L,(v) we need to worry about what 
sets the scale for uo . As we have discussed above , this has to be the rate of 
decay of fluctuations r(£). In the finite geometry that we are discussing now 
, the scale for decay of fluctuations will depend on L as well . In discussing 
the correction depicted in Eqn(5) , it is obvious that this fine point need 
not be discussed as this correction is already O(^) . For He (superfluid 
transition) , there is in someways an additional simplifying feature . For 
the order parameter decay rate the non-linear effect of fluctuations becomes 
significant , only very close to the critical point and for all practical purposes 
, the relaxation rate can be taken to be at its non critical background value . 

The complex order parameter field ipi(x)i = 1,2 will be governed by the 
Langevin equation 

5F 

1>i = -r — + Ni 6 

where 

F = jd D x[^ + \wy + ^r} (7) 
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and N is a Gaussian white noise . For reasons stated above we choose to 
drop the reversible term (the Josephson equation for the phase of the order 
parameter) . The parameter m 2 is proportional to T — T A , where T A is 
the bulk transition temperature . The system is confined in one of the D 
directions . We call that the z-direction . It is convenient to work with 
the fourier transform in D — 1 directions and the fourier series (Dirichlet 
boundary conditions at z = and z = L suppressing the fluctuations) in the 
z direction . The expansion of the time-dependent order parameter field is 

Ur,t)=J2Mn,K,t)exp* K - R sm(—) (8) 

n L 

The equation of motion for ipi(n, K, t) is 

Tl 2 7Y 2 

^(n, K, t) = -r (m 2 + K 2 + -^-)^(n, K, t) + N t + 0(^ 3 ) (9) 

In what follows , we will assume that all static correlations have been ac- 
counted for and m 2 = £~ 2 . The specific heat is obtained as the response 
function corresponding to the time dependent correlation function 

D(£, L, t 12 ) = JII d Zl dz 2 d D R 12 < V> 2 (i?i, Zl , ti)^ 2 (i? 2 , Z2, t 2 ) > (10) 

with D(£,L,u) = 2 ImC ^' L '^ according to fluctuation dissipation theorem , 
straightforward algebra leads to (two term accuracy , L — > oo limit and the 
first correction) 

d D p 1 1 



C(t,L,u) 



(2tt) d (p 2 + m 2 ) (-^ + p 2 + m 2 ) 

1 r d D ~ x p 1 1 
2L J (2tt) d - 1 (p 2 + m 2 ) (-^ + p 2 + m 2 



(11) 



We work to logarithmic accuracy and hence evaluate the integrals at D = 4 
(proper exponentiation can be undertaken by working to two loop order , 
the details of which will be published elsewhere) to get the functions f(Q) 
and a(fi) introduced in Eqs(3) and (5) . Note that since we are taking the 
logarithmic divergence for the bulk specific heat , the C(£) in Eqns(4) and 
(5) reduces the constant C of Eqn(3) . The function f(fl) and a(Q) are 

/(fi) = I(-L-l)ln(l-^) (12) 
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a(n) = ^-^[VT=m-l] (13) 



leading to 



C(£,L,u;) = C {ln— -^ln(l + fi 2 )--^-tan^ 1 (fi)+i[itan' 1 (fi) 

1 , ^9 Xn 7T ^ 9 .i . ,tan _1 n. in 

■ ln(l + Q 2 )} - — — (1 + fi 2 ) 4 sm(^— ) - 



4fi v /J mLQ x ' y 2 ' mm 

r , H 1 / "t . -, 

[ ( 1 + ^ )3C0S (__)_1] } 

= C7* + »Cj (14) 

where Cr and Cj are the real and imaginary parts of the specific heat. Con- 
sidering the zero frequency limit we see that Cr = Co[mA/m — n/2mL] 
leading to a = n/2 in Eq. (4) which is to be compared with a ~ 1.4 obtained 
in [11]. 

We now return to Eqn(l) , to find the attenuation and dispersion . The 
attenuation per wavelength is |£ = ^ffi^-, which leads to the fre- 
quency attenuation {uj ^> 2r m 2 ) as 



2vr Mo [ln(^) - v^^] 2 + [i _ 2v ^(^ 



(15) 



This is the saturation attenuation per wavelength , which does not change 
as the temperature is lowered further , where u>o/2n = 30GHz,T = 1.2 x 
10- 4 cm 2 sec-\ ui/uo = 8/3 x 10" 2 . 

For the plate separation of 211CL4 of mehta and Gasparini , the reduc- 
tion in the attenuation due to the quenching of fluctuations is about 18% at 
10MHz and increases to 45% at 2.5MHz . This is a large effect compared 
to the 4% surface effects that show up in the static measurements . For the 
corresponding measurement of thermal conductivity near the superfluid tran- 
sition Kahn and Ahlers [9] found that the deviation from the bulk is about 
7% when the correlation length £ equals the confining length L (in their case 
the radius of the pore) . The surface effect for the ultrasonic measurement 
can easily amount to 30% which makes this an attractive system for a con- 
frontation between theory and experiment . The FSE on the dispersion can 
be obtained from the real part of Eqn(l). 
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We note that the above is a one loop calculation in the critical region . 
The lack of crossover to the background in our treatment of the specific heat 
implies that we can consider frequency uj which are much smaller than the 
cut-off frequency uj . This is a restriction on the validity of the dashed curve 
shown in Fig(l) . The solid curve in addition is restricted to confining lengths 
which are not too small i.e. L > (2r /c<;)2 and in this regime the accuracy 
of the calculation is restricted by the loop order . This is not too severe 
a restriction as an accuracy of 0(e) which our calculation entails,becomes 
an accuracy of 0(f) when the combinatorial factors are included . Thus, 
in the above mentioned ranges of the parameters uj and L , the dashed 
curve in Fig(l) should be an accurate prediction . It should be noted that 
contrary to the static specific heat or the thermal conductivity , the sound 
properties can only be probed in a real experiment . The final issue , then 
, is whether the effect can be observed in real a experiment . The critical 
ultrasonics near the superfluid transition has been studied more than two 
decades ago . The most accurate data of that period lies in the 0.5MHz — 
5MHz range . In this region the scatter in the data is about 15%. This is 
somewhat better than the borderline for detecting the suppression reported 
here . Considering the fact that developments in the experimental field would 
enable more accurate measurements at the present time , we believe this effect 
should be experimentally accessible . 

The other sensitive part of an ultrasonic measurement is the low fre- 
quency end {uj <C 2r m 2 ) , where for the bulk substance the attenuation per 
wavelength is proportional to C 2 R Vl/A . The relative correction for the FSE 
is 1 — 2^ , once again a larger effect than can be obtained in statics. For an 
easily realizable situation of ml ~ 8 this gives a 20% reduction in the attenu- 
ation . The whole course of the attenuation function with its dependence on 
uj and L is straightforward to obtain and will be exhibited elsewhere . Here 
we have reported the salient feature , which carry the most experimentally 
accessible signatures . We hope this will stimulate experimental activity in 
the field. 
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Figure caption 

Fig.l. Saturation attenuation is plotted against frequency.The dashed 
curve shows the bulk (L — > oo) result whereas the solid curve shows the 
surface effect. 
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